Bound i7-dibaryon in Flavor SU(3) Limit of Lattice QCD 
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The flavor-singlet _ff-dibaryon, which has strangeness —2 and baryon number 2, is studied by the 
approach recently developed for the baryon-baryon interactions in lattice QCD. The flavor-singlet 
central potential is derived from the spatial and imaginary-time dependence of the Nambu-Bethe- 
Salpeter wave function measured in Nf = 3 full QCD simulations with the lattice size of L ~ 2, 3, 4 
fm. The potential is found to be insensitive to the volume, and it leads to a bound H-dibaryon with 
the binding energy of 30-40 MeV for the pseudo-scalar meson mass of 673-1015 MeV. 

PACS numbers: 12.38. Gc, 13.75.Ev, 14.20.Pt 



Search for dibaryons is one of the most challenging 
theoretical and experimental problems in the physics 
of strong interaction and quantum chromodynamics 
(QCD). In the non-strange sector, only one dibaryon, the 
deuteron, is known experimentally. In the strange sector, 
on the other hand, it is still unclear whether there are 
bound dibaryons or dibaryon resonances. Among others, 
the flavor-singlet state (uuddss), the i7-dibaryon, has 
been suggested to be the most promising candidate [H. 
The H may also be a doorway to strange matter and to 
exotic hyper-nuclei Although deeply bound H with 
the binding energy Bh > 7 MeV from the AA threshold 
has been ruled out by the discovery of the double A nu- 
clei, A^He f3|, there still remains a possibility of a shallow 
bound state or a resonance in this channel 

While several lattice calculations on H have been re- 
ported as reviewed in Q (see also recent works [B-[l|), 
there is a serious problem in studying dibaryons on the 
lattice: To accommodate two baryons inside the lattice 
volume, the spatial lattice size L should be large enough. 
Once L becomes large, however, energy levels of two 
baryons become dense, so that quite a large imaginary- 
time t is required to make clear isolation of the ground 
state from the excited states. All the previous works on 
dibaryons more or less face this issue (see also [§|). 

The purpose of this Letter is to shed a new light on 
the i/-dibaryon by extending the lattice approach re- 
cently proposed by the present authors Our start- 



ing point is the baryon-baryon potential obtained from 
the Nambu-Bethe-Salpeter (NBS) amplitude measured 
on the lattice [13] ■ Such a potential together with the 
NBS amplitude can be shown to satisfy the Schrodinger 
type equation and to reproduce the correct phase shifts 
at low energies. It was found on the lattice in the flavor 
SU(3) limit [7| that, while the celebrated repulsive core of 
the potential appears in the nucleon-nucleon(NN) chan- 
nels, the "attractive core" emerges in the _ff-dibaryon 
channel. These features at the short range part of the 
potential are essentially dictated by the Pauli exclusion 
principle in the quark level: Six-quarks residing at the 
same spatial point is partially forbidden by the quark 
Pauli effect in the NN channels, which belong to the fla- 
vor 27-plet or 10*-plet, while the flavor-singlet six-quarks 
do not suffer from the Pauli effect 11 1 (see also |12{). 



The approach based on the baryon-baryon potential 
has several advantages. In particular it can be used not 
only to reduce the finite volume artifact but also to avoid 
the problem of contaminations from excited states, as will 
be explained later. In this Letter, to capture essential fea- 
tures of the H-dibaryon without being disturbed by the 
quark mass differences, we consider the flavor SU(3) limit 
where all u, d, and s quarks have a common finite mass. 
This allows us to extract baryon-baryon potentials for ir- 
reducible flavor multiplets and to make the comparison 
among different flavor channels in a transparent manner. 
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We start with the NBS wave function flOl defined by 

0„(f) = (O|(i3i?)(")(f,O)|Ty„;a), (1) 

where the state vector |W^„; a) is a QCD eigenstate with 
the baryon number 2 (6 quark state) and energy Wn 
in the flavor a-plet. {BB)'^°'\r,t) = Y.r,j,sC\f B,{x + 
r, t)Bj(x, t) is a two-baryon operator with a relative dis- 
tance f in a-plet with Bi being a one-baryon composite 
field operator in the flavor octet. The relation between 
two-baryon operators in the flavor basis and baryon basis 
are given by the SU(3) Clebsch-Gordan coefficients. 

In the lattice QCD simulations, the above NBS wave 
functions is extracted from the four point function as 
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Here (BB) (to) is a wall source operator at time to to 
create two-baryon states in a-plet, while (iJB)'^"^ (r, t) 
is the sink operator at time t to annihilate the two- 
baryon states. Even if we choose t — to moderately 
large so that the inelastic scatterings (e.g. the scatter- 
ing with excited baryons and the scattering with me- 
son production) do not contribute to G4, there still re- 
main elastic scattering states with low energy excitations 
due to the relative motion of the baryons. For exam- 
ple, with the baryon mass M ~ 2 GeV in a finite box 
of i = 4 fm, the non-interacting two-baryon system 
has Wi-Wo ^ {2tt / Lf / {2y) ~ 50 MeV, with the re- 
duced mass ^ — M/2. This requires t — to > 10 fm to 
achieve 1/10 suppression of the first excited state 0i(r) 
in Gi{f,t — to). It is beyond most of the previous and 
current lattice simulations. 

Our potential approach avoids the above problem 
in the following way: The two-body potential in low 
energy QCD dictates all the elastic scattering states 
0„(r, t) = (/)„(r)e~^'^"~^*^)' simultaneously through the 
Schrodinger equation in the Euclidean space-time [lo| . 
With the non-relativistic approximation for W„ , it reads 

HoMr,t)+ I d'r'Uir.T^Mi' ,t) ^ ^^6r,{T,t), (3) 



dt 

where Ho = — V^/(2/i) and [/ is a non-local and energy- 
independent potential. Since the above equation is lin- 
ear in (f)n, the linear combination such as (t>{'r,t) = 
J2n^n4>nif,t) = G4(r,t)/e-2A" also satisfies Eq.®. 
We note that the derivative expansion of U in terms of 
its non-locality leads to U{r,r') = [Vc(f) -I- VT{r)Si2 + 
VLs{r)L-S+-- ■ )5{r~r') 0, where Vc, Vt and Vls are 
the central, tensor and spin-orbit potentials, respectively, 
and dots stands for terms including power of V. It was 
shown in [3] that the leading order potentials without 
V dominate the potential at low energies. Thus, the rel- 
evant term in the spin-singlet channel, Vc^ is obtained 



TABLE I: Summary of lattice parameters and hadron masses. 
The uncertainty of a 18] is not reflected in hadron masses. 



a [fm] 


L [fm] 


^uds 


mp, [MeV] 


ruB [MeV] 


^cfg 






0.13710 


1015.0(6) 


2030(2) 


360 


0.121(2) 


3.87 


0.13760 


836.5(5) 


1748(1) 


480 






0.13800 


672.9(7) 


1485(2) 


240 
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In this way, one can extract the baryon-baryon potential 
without identifying each elastic states cj)n{r,t) as long 
as t — to is so chosen that the inelastic scatterings are 
suppressed. Once we obtain the volume independent Vc, 
binding energies and scattering phase shifts in the infinite 
volume are obtained by solving the Schrodinger e quat ion. 
In contrast to the conventional Liischer's method [IJ], we 
do not calculate the energy shift of two hadrons at finite 
L to access the observables at L — >■ 00. Further theo- 
retical details of this method will be given in a separate 
publication 15 1. 

Let us now consider the interaction between flavor- 
octet baryons in the flavor SU(3) limit, for which two 
baryon states with a given angular momentum are la- 
beled by the irreducible flavor multiplets as 8 ® 8 = 

(27 © 8s © l)syminetric ® (10* ffi 10 © 8a)anti-syinmotric- 

Here "symmetric" and "anti-symmetric" stand for the 
symmetry under the flavor exchange of two baryons. 
For the system in the orbital S-wave, the Pauli prin- 
ciple between two baryons imposes 27, 8^ and 1 to 
be spin singlet (^S'o) while 10*, 10 and 8^ to be spin 
triplet (^Si). Since different multiplets are independent 
in the flavor SU(3) limit, one can define the correspond- 
ing potentials as ^(^'^^(r), V^^'\r), y(i)(r) for ^^o and 
Y{W)(^r), y(i°)(r), l/(«'')(r) for 3S'i. Hereafter, we fo- 
cus on the flavor-singlet channel with 




(5) 



where A, E, N and S are the standard baryon operators 
with Lorentz structure, [qiC^5)q]q 

In our dynamical lattice QCD simulations, we employ 
the renormalization group improved Iwasaki gauge action 
and the non-perturbatively 0{a) improved Wilson quark 
action. For 16^ x 32 lattice, we use the configuration set 
generated by CP-PACS and JLQCD Collaborations fll] 
aX (3 — 1.83. In addition, we generate gauge configura- 
tions with the same /3 for 24'^ x 32 and 32'^ x 32 lattices, 
using the DDHMC/PHMC code [l3|. Quark propaga- 
tors are calculated for the spatial wall source at to with 
the Dirichlet boundary condition in the temporal direc- 
tion. The sink operator is projected to the represen- 
tation of the cubic group, so that the NBS wave function 
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FIG. 1: Flavor 27-plet potential V^^ (r) obtained for lattice 
sizes L = 1.94, 2.90, 3.87 fm at trips = 1015 MeV and {t - 
to)/a = 10. 



is dominated by the S-wave component. For the time 
derivative, we adopt the symmetric difference on the lat- 
tice. Lattice parameters such as lattice spacing a, the 
hopping parameter K„ds, the number of configurations 
iVcfg, together with the pseudo-scalar meson mass mpg 
and the octet baryon mass ms are summarized in Table 
Ufor 323 X 32 lattice. 

To check the qualitative consistency with previous 
works, we show in FiglT] the central potential in the 27- 

(27) 

plet channel (^) obtained in three different lattice 
volumes with L ^ 1.94, 2.90, 3.87 fm at rrips = 1015 
MeV and {t — to)/a = 10. This is the case correspond- 
ing to the NN potential in the ^Sq channel. Compared 
with statistical errors, the L dependence is found to be 
negligible. The t dependence is also small as long as 
{t — to) /a > 9. Note that we do not need overall shift 
of the potential: it approaches zero automatically as r 
increases. The figure shows a repulsive core at short dis- 
tance surrounded by an attractive well at medium and 
long distances, which is qualitatively consistent with our 
previous results in quenched and full QCD simulations 
reviewed in [lij . 

Shown in Figl2l[a) and Figl2][b) are the volume depen- 
dence and the quark mass dependence of the central po- 
tential in the flavor-singlet channel V^'^\r), respectively. 
In both figures, we take {t — to)/a = 10 and have checked 
that the potentials do not have appreciable change with 
respect to the choice of t. We find that the flavor-singlet 
potential has an "attractive core" and its range is well lo- 
calized in space. Because of the latter property, we find 
no significant volume dependence of the potential within 
the statistical errors as seen in FiglSfa). We find that the 
long range part of the attraction tends to increase as the 
quark mass decreases [FigHIb)]. 

We fit the resultant potential by the following analytic 
function composed of an attractive Gaussian core plus 
a long range (Yukawa)^ attraction: V{r) — feie^^^'' + 
63(1 — 6"^*"'^) (e~''^''/r) . With the five parameters, bi - 
&5 , we can fit the function to the lattice results reasonably 
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FIG. 2: Flavor-singlet potential V^\r) at {t - to) /a = 10. 

(a) Results for L = 1.94, 2.90, 3.87 fm at mps = 1015 MeV. 

(b) Results for L = 3.87 fm at rrips = 1015, 837, 673 MeV. 



weU with xVdof ~ 1. The fitted result for L = 3.87 fm 
is shown by the dashed line in Fig 121 a). 

Finally, using the potential fitted by the function, we 
solve the Schrodinger equation in the infinite volume and 
obtain the energies and the wave functions for the present 
quark masses in the flavor SU(3) limit. It turns out that, 
in each quark mass, there is only one bound state with 
the binding energy of 30-40 MeV. In FiglHlJa), the energy 
and the root-mean-square (rms) distance of the bound 
state are plotted in the case of {t ~ to)/a = 9, 10, 11 at 
TOps = 673 MeV and L ~ 3.87 fm, where errors are esti- 
mated by the jackknife method. Although the statistical 
error increases as t increases, we observe small changes 
of central values, which will be included as the system- 
atic errors in our flnal results. FigjHJb) shows the energy 
and the rms distance of the bound state at each quark 
mass obtained from the potential with L — 3.87 fm and 
{t — to)/o. — 10. Despite that the potential has quark 
mass dependence, the resultant binding energies of the 
i7-dibaryon are insensitive in the present range of the 
quark masses. This is due to the fact that the increase of 
the attraction toward the lighter quark mass is partially 
compensated by the increase of the kinetic energy for 
the lighter baryon mass. It is noted that there appears 
no bound state for the potential of the 27-plet channel in 
the present range of the quark masses. 

The final results of the binding energy in the SU(3) 
limit Bh and the rms distance -y/ (r^) are given below. 
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m„ =1015 [MeV] — e- 
mps = 837 [MeV] 
m,s = 673 [MeV] 
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FIG. 3: Bound state energy Eq = —Bh and the rms distance 
\/Jr^ of the //-dibaryon obtained from the potential at L = 
3.87 fm. (a) Imaginary-time dependence at mps — 673 MeV. 
(b) Quark mass dependence at (t — to) /a — 10. 



with statistical errors (first) and systematic errors from 
the i-dependence(second). 



m„s ^ 1015 MeV : Bh = 



mps = 837 MeV : Bh 



mps = 673 MeV : Bh 




32.9(4.5)(6.6) MeV 
0.823(33) (40) fm 
37.4(4.4)(7.3) MeV 
0.855(29)(61) fm 
35.6(7.4)(4.0) MeV 
1.011(63)(68) fm 



A less than 1% error from the choice for the fit function 
is not included here. 

Since Bh has the weak quark mass dependence, one 
may assume a similar binding energy is realized even with 
the realistic SU(3) breaking, where Bh is interpreted as 
the binding energy from the average mass of two octet 
baryons in the S = —2 and / = channel. Considering 
that the difference between this average and 2mA is about 
the same amount to Bh, the _ff-dibaryon may appear 
as a weakly bound state or a resonant state near the 
A A threshold, as mentioned in [3]. To make a definite 
conclusion on this point, however, we need (2 -1-1) -flavor 
lattice QCD simulations with the AA — A^S — SS coupled 
channel analysis as well as a careful study on the non- 
locality of the potential. The extension of the method 
outlined in this Letter to this direction is in progress [l^] . 
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